arXiv:1506.01484v4 [quant-ph] 19 Oct 2016 


Evaluation of entanglement measures by a single observable 

Chengjie Zhang/’^’Q Sixia Yu,“’^ Qing Chen,“’"^ Haidong Yuan,^ and C.H. Oh“’® 

'College of Physics, Optoelectronics and Energy, Soochow University, Suzhou, 215006, China 
^ Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore 
'Hefei National Laboratory for Physical Sciences at Microscale and Department of Modern Physics, 
University of Science and Technology of China, Hefei, Anhui 230026, China 
‘'Department of Physics, Yunnan University, Kunming, 650091, China 
'Department of Mechanical and Automation Engineering, The Chinese University of Hong Kong, Hong Kong 
'Physics Department, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore 

(Dated: October 20, 2016) 

We present observable lower bounds for several bipartite entanglement measures including entanglement 
of formation, geometric measure of entanglement, concurrence, convex-roof extended negativity, and G- 
concurrence. The lower bounds facilitate estimates of these entanglement measures for arbitrary finite¬ 
dimensional bipartite states. Moreover, these lower bounds can be calculated analytically from the expectation 
value of a single observable. Based on our results, we use several real experimental measurement data to get 
lower bounds of entanglement measures for these experimentally realized states. In addition, we also study the 
relations between entanglement measures. 


I. INTRODUCTION 

Quantum entanglement is widely recognized as a valuable 
resource in quantum information processing. However, it is 
far from simple to fully determine entanglement. Therefore, 
the characterization and quantification of entanglement be¬ 
come fundamental problems in quantum information theory. 
Lots of entanglement measures have been proposed, such as 
entanglement of formation, geometric measure of entangle¬ 
ment, concurrence, and convex-roof extended negativity. 

Consider a finite dimensional bipartite system, one is sub¬ 
system A and the other one is subsystem B. The entanglement 
of formation (EOF) is the first entanglement measure built by 
the convex roof construction iSli. For a pure state |i/r), it 
is defined by Ef{\ij/)) - S{qa), where S{q) - -Tr(plog 2 p) 
stands for the von Neumann entropy and qa - TrB(|i/r)(t/r|) 
is the reduced density matrix of subsystem A. For a mixed 
state Q, the EOF is defined by the convex roof, Ef(g) = 
infjp; !,/,;)) 2; p/Ff (|i/r,)) for all possible ensemble realizations 
Q - where p, > 0 and Xi Pi - 1- The infi- 

mum represents the minimal possible average entanglement 
over all pure state decompositions of g. The geometric mea¬ 
sure of entanglement (GMK) is another kind of convex-roof 
entanglement measures ll^01- For an arbitrary bipartite pure 
state \i//) - Ua ® UbYji with s/jui being its Schmidt 

coefficients, the GME is defined by Fed'/')) = 1 “ rnax{p,). 
Similarly, the GME is extended to mixed states by the con¬ 
vex roof. The concurrence was firstly introduced for two- 
qubit states by Ref. H, and based on it Wootters and co¬ 
workers derived computable formulae for concurrence and 
EOF in the two-qubit case lHll]. After that, Refs. 0- 
lit] extended it to bipartite higher-dimensional systems, i.e., 

= .^2(T^-TYp^ for pure states, and its convex roof for 
mixed states. The convex roof extended negativity (CREN) is 
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generalized from the negativity, which is strongly related with 
the partial transpose lUOll . For a bipartite state g, its negativ¬ 
ity reads Nig) - ||p^®|| - 1 (for simplicity we ignore the 
coefficient 1 /2 ), where || - || stands for the trace norm and Tb 
is partial transpose with respect to subsystem B. The positive 
negativity is a necessary and sufficient condition of entangle¬ 
ment for pure states, 2x2 and 2x3 mixed states, but only a 
sufficient condition for higher-dimensional mixed states lUdl . 
To overcome this drawback, Lee et al. proposed the CREN 
id. For a pure state CREN is defined by the nega¬ 
tivity Ni\if/)) - N{\4i)) - |||tA)(iAr®ll - L For mixed states, 
CREN is defined by the convex roof. Last but not least, the 
G-concurrence can be defined by the determinant of the re¬ 
duced density matrix: (2q(|iA)) = m(detpA)*^'" for an m ® « 
{m < n) pure state \tft) ol [o [f^ . and for mixed states, the 
G-concurrence is defined by the convex roof as well d. 

Although many entanglement measures have been pro¬ 
posed fl Ul^ . there are only a few explicit expressions of 
these measures for two-qubit states and some ^cial kinds 
of higher-dimensional mixed states flla [lZ[l6l - [l^ . Further¬ 
more, for a general state it is proved that computing many en¬ 
tanglement measures including the entanglement of formation 
is NP-hard |d.|2^, which implies that we could only derive 
bounds on (rather than compute exact values of) these entan¬ 
glement measures. Thus, in order to evaluate entanglement 
measures, lower and upper bounds of entanglement measures 
for general higher-dimensional states and multipartite states 
have been proposed ll2TI - [^ . Besides, the lower bounds of 
entanglement measures can server as a valuable tool for opti¬ 
mal control lH^. However, if the proposed bound cannot be 
directly measured in experiments, quantum state tomography 
has to be performed which leads to rapidly growing experi¬ 
mental resources as system size increases. Therefore, experi¬ 
mentally observable lower and upper bounds of entanglement 
measures attract much interest recently 

In this paper, we shall propose observable lower bounds 
for EOF, GME, concurrence, CREN, and G-concurrence in 
finite-dimensional bipartite systems. These lower bounds can 
be obtained from the expectation value of a single observable. 
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Based on our results, we will present several examples using 
real experimental measurement data. Furthermore, the rela¬ 
tions between entanglement measures will be studied. 


II. ENTANGLEMENT OF FORMATION 


For simplicity, we use the denotation co(g). Here co(.g) de¬ 
notes the convex hull of the function g, which is the largest 
convex function that is bounded above by the given function 
g. The denotation has been used to get explicit expressions 
and bounds for the EOF IHraisolM. 

Theorem 1 .— For any m® n {m < n) quantum state g, its 
entanglement of formation Ef(g) satisfies 

Epig) > co[f?(A)], ( 1 ) 


Schmidt coefficients in decreasing order, and we have used 



m sitn 


{4>\Ui®U2\ilrj){ilrj\ul®ul\<l>) 

> max-, ( 4 ) 

Ui.Ui s\m 


where the second inequality holds since the theorem shown 
in Ref. lHH], and the detailed proof has been given in the ap¬ 
pendix. Therefore, 




max 

t/l,C/2 


S\m 


S\m 


where A = max{((/i|p| 0 )/(sim), 1 /m}, R(A) = H2[y(A)] H- [1 - 
■y(A)] log2(m - 1 ), 7(A) = [ VX + y/(m - 1)(1 - A)]^/m, with 
H2ix) = -xlog2X - (1 - A:)log2(l - x) being the standard 
binary entropy function, and \<p} -Va^Vb TiT=i VX|b) being 
an arbitrary pure entangled state in m®n system (where { yisv) 
are its Schmidt coefficients in decreasing order). The convex 
hull of R(A) is 


co[R(A)] 


H2[y{^)] + [1 - r(A)] log2('« ■ 

Ag [1 4 ^] 


1 ), 


2li£||^(A-l) + log2m, 

Ae 

V L , 11 - ’ J 


( 2 ) 


Together with T-' = iYIJLi one can get 

YjjPj'^^ ^ max{{ 0 |p|(^)/(iim), 1 /m) = A. Since co[R(/l)] is 
a monotonously increasing function and Pj'^^ ^ one 
has co[R(Yj Pj'^^)] S co[R(A)]. 

As introduced above, co(g) is the largest convex function 
that is bounded above by the given function g. From the ex¬ 
pression of R(A), the explicit expression of co[R(A)] is as Eq. 
(I 2 ]i lU^ 1 ^ pll] . Actually, the same function has been got¬ 
ten for the EOE of isotropic states S and the lower bound 
of EOE based on partial transpose and realignment criteria in 
Refs. miH. ■ 


Proof- We first find the minimal admissible H{fj) = 
-2™iA'/log2A'/ = -Tr(plog2p) = 5 (p) for a given d = 
{Y 7 =i yfPif 'where /r, are eigenvalues of g and fl is the 
Schmidt vector {//i,/i2, • • • ,Pm]- Consider the following func¬ 
tion, 

R{A) = mm|//(/l)d = |’ 


III. GEOMETRIC MEASURE OF ENTANGLEMENT 

Similar to entanglement of formation, we can also find an 
observable lower bound for the GME. 

Theorem 2 .— Eor any m ® n (m < n) quantum state g, its 
geometric measure of entanglement Ecig) satisfies 

Ec(g) > co[e(A)], (6) 


As shown in Ref. Cl, the minimal Hip) versus A corresponds 
to p in the form {f, (1 - t)/{m - 1), • • • , (1 - t)/(m - 1)) for 
t e [ 1 /m, 1 ]. Therefore, one can get R(A) - H2[y{A)] -H [1 - 
7(d)] log2(m - 1) with 7(d) = [ Vd -H p{m - 1)(1 - d)]^/m. 

Suppose that we have already found an optimal decompo¬ 
sition YjPM>pi'i'i\ for Q to achieve the infimum of Ef(g), 
then ^^(p) = by definition. Since co[R(d)] 

is a monotonously increasing convex function and satisfies 
co[R(d)] < R{A) < H(p) for a given d, one thus has 

Efig) = = 

J J 

^ ^Pfo[R(A^)] > co[R(^p2d^)] > co[R(A)], 

j j 

where \ipj} ^ Ua ® UbYZi with {ypP being its 


where A = max{((/i|p| 0 )/(sim), 1 /m), 2 (A) = 1 - 7(A) with 
7(A) = [ VA + p{m — 1)(1 - A)]^/m, and 

co[e(A)] = e(A). (7) 

Proof .— We first find the minimal admissible G{p) s 1 - 
gLmax for a given d = (2” j yfp'ifim, where g.max = max{//,). 
Consider the following function. 


2(d) = min 

P 




( 8 ) 


As shown in the appendix, the minimal G{p) versus d corre¬ 
sponds to p in the form (f, (1 - t)l{m - 1), • • • , (1 - t)l{m — 1)) 
for t e [ 1 /m, 1 ). Therefore, 2 ('f) = 1 “ t(^) holds, with 
7(d) = [ Vd + p(m — 1)(1 - d)]^/m. 

Similar to the proof of Theorem 1 , suppose that we have 
already found an optimal decomposition Yj Pj\'l^j)('l^j\ for p 
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to achieve the infimum of Ec(q), then Ea(g) — YjJ PjEc{\4'j)) 
by definition. Since co[2(/l)] is a monotonously increasing 
convex function and satisfies co[2(/l)] < Q(A) < G(jj) for a 
given A, one thus has 

Ec(e) = ^PjEcQtf'j)) ^^PjGip-') 

J J 

> YjPjCo[Q(Aj)] > co[Q(J]pjAj)] > co[e(A)], 

J J 

where, similar to the proof of Theorem 1, we have used 
I^jPjAj > A. 

From the definition of co, one can see co[2(A)] = 2(A), 
since Q{A) is a convex function. The same function has been 
used for the GME of isotropic states H. ■ 

IV. CONCURRENCE 

In the following, we shall also present an observable lower 
bound for the concurrence. 

Theorem 3. — For any m®n{m < n) quantum state g, its 
concurrence C{g) satisfies 

C(g) > co[P(A)], (9) 

where P(A) = -^2[ 1 - yC A)] [my( A) + m-2\I (m - 1), 

with 7 (A) = [VA + V(™ “ 1)(1 “ A)]^/m, A = 

max{{<p\g\(l))/(sim), 1 /m), and 


co[P(A)] = J^^^(A-i). (10) 

Proof. — Similarly, we first find the minimal admissible 

E{p) = for a given A = (ZT=i Con¬ 

sider the following function, 

P{A) = n:un|L(/l)d = j- (H) 

As shown in the appendix, the minimal L{fj) versus A cor¬ 
responds to yu still in the form {f, (1 - t)l{m - 1), ••• ,(1 - 
f)/(m - 1)} for t e [1/m, 1]. Therefore, P(A) - 
-y^2{T^^7p)]|m7p)Tm^^'2i7(m^n3 holds, with y(A) - 
[VI+ i(m-l)(l-A)f/m. 

Similar to the proofs of Theorem 1 and 2, suppose that we 
have already found an optimal decomposition YjJ Pj\^j)(^j\ 
for g to achieve the infimum of C(g), then C{g) - 
Yjj PjC{\tl/j)) by definition. Since co[P(/l)] is a monotonously 
increasing convex function and satisfies co[P(/l)] < L(jj) for a 
given A, one thus has 

C(g) - ^/7;C(|<A,» = 2p;T(/?^) 

J J 

> Yj PjCo[P{AJ)] > co[P(^ pjA^)] > co[P(A)], 
j J 

where 2; PjA^ > A has been used again. 


From the definition of co, one can get Eq. Coll, since P{A) 
is a monotonously increasing concave function as shown in 
the appendix. The same function has been gotten for the 
concurrence of isotropic states dl, the lower bound of con¬ 
currence based on partial transpose and realignment criteria 
in Ref. and the special case of Eq. (l9]l with \4>) - 

1 / yfinYUli |h) being the maximally entangled state 11 ^D^ . 


V. CONVEX-ROOF EXTENDED NEGATIVITY 

Similarly, an observable lower bound of CREN has been 
presented as follows. 

Theorem 4 .— Eor any m ® n (m < n) quantum state g, its 
convex-roof extended negativity M(g) satisfies 

N{g)>mA—\, (12) 

where A = max{(</i|p|0)/(iim), 1/m). 

Proof .— Similar to the proofs of the theorems above, sup¬ 
pose that we have already found an optimal decomposition 
pj\\l/j){\l/j\ for g to achieve the infimum of N{g), then 
El{g) - Yj PjEfilfj)) by definition. It is worth noticing that, 

for an arbitrary pure state = Ua^UbYi with fri 

being its Schmidt coefficients in decreasing order, we have 

= (Zti fp{f - 1- Thus, 

rn —— 

Me) = 

j j '■=! 

= m pjA-' — 1 > mA — 1 , 
j 

where Yj Pj'^^ ^ ^ has been proved in Theorem 1 . Actually, 
a similar function has been gotten for the CREN of isotropic 
states ifl^ . ■ 


VI. G-CONCURRENCE 


Last but not least, one can present an observable lower 
bound of G-concurrence as follows. 

Theorem 5 .— Eor any m® n {m < n) quantum state g, its 
G-concurrence Ccig) satisfies 

Ccio) > co[MA)], (13) 

where fi:(A) = m[ 7 (A)j 6 (A)"'-']h'", with 7 (A) = [ Va- 
^J{m-V){l-A)flm,l3(A) = [ VX -H V(1 - A)/(m - V)flm, 
A - m?ix{{(j)\g\(p)l{s\m), 1/m), and 

co[fi:(A)] = maxjl - m(l - A),0). (14) 

Proof .— Similarly, we first find the minimal admissible 
Sift) = for a given A = (Yti 'JPifim. Con¬ 

sider the following function. 


K{A) - rmn-[5(ju) 



(15) 
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As shown in Ref. if^ and the appendix, the min¬ 
imal S (ft) versus A corresponds to fl in the form 
{f,(l - t)l(m - - t)l(m - 1 )) for t e 

[0, l/m]. Therefore, K(A) = m[y(A)l}(A)"'-^YI”' holds, with 
y(A) - [yfl - V('« - 1)(1 - and ft{A) - [ VI ^- 

Similar to the proofs of above theorems, suppose 
that we have already found an optimal decomposition 
Yjj Pj\>fj){^j\ for to achieve the infimum of Ccig), then 
Cc(g) - Zij PjCcil^j)) by definition. Since co[A'(/f)] is 
a monotonously increasing convex function and satisfies 
co[fir(/l)] < S(ft) for a given A, one thus has 

Ccig) = = 

J J 

^ ^Pfo[K(A^)] > co[K{^pjAJ)] > co[/:(A)], 

J j 

where 2 ; Pj'^^ ^ A. has been used again. 

From the definition of co, one can get Eq. (fl4li . since 
K{A) is a monotonously increasing concave function in [{m - 
1)/ot, 1] as shown in Ref. ifs^ . The same function has been 
gotten for the G-concurrence of axisymmetric states ifs^ . and 
the lower bound of G-concurrence for arbitrary states which 
is the special case of Eq. (fT3l l with \^) - 1 / ^|m YiT=i I”) being 
the maximally entangled state if^ . ■ 

Remark 1 .— It is worth noticing that when we choose a spe¬ 
cial case for | 0 ), i.e. \(p) - |iA^> where \ip^) - 1 / \ii} 

is the maximally entangled state, then becomes to l/m 
and A = max{{if^\g\if^),l/m}. Eor this special case, our 
results shown in all the above theorems can be proved in a 
different manner: for an arbitrary m ® m state g, one can 
project it onto the isotropic states by the twirling opera¬ 
tion, i.e. p'™ = J dU(U ® U*)g(U ® t/*)V which is a lo¬ 
cal operations and classical communication (LOCC) opera¬ 
tion and therefore cannot increase entanglement. Moreover, 
we have (iA^|p|iA^) = since it is invariant un¬ 

der the twirling operation. The entanglement measures of the 
isotropic states can be expressed as a function of 
1141 [la [la [l^ . Thus, one can get lower bounds of entangle¬ 
ment measures for g from the entanglement measures of the 
isotropic states, this idea to get a lower bound has been known 
since the earliest paper |[ll]. However, this alternative proof for 
10 ) = is not valid for a general \<p). 

Remark 2 .— Erom Theorem 1 to Theorem 5, all the lower 
bounds proposed above are the functions of A, which only 
depends on the expectation value of a single observable | 0 ){ 0 |. 
Therefore, it will be much easier to evaluate than tomography 
in experiments. 


VII. EXAMPLES AND EXPERIMENTAL 
MEASUREMENTS 


The first example is a real experimental state shown in 
Ref. S. Tonolini et al. experimentally realized a 
high-dimensional two-photon entangled state, with dimen¬ 
sion of each photon being equal to d - 17. They recon¬ 
structed the density matrix p„pj of the experimental state. 


and found the fidelity with the maximally entangled pure 
state being Tx ^ ^,/g^ = 0.831. Therefore, 

our parameter for this experimental state should be Ai = 
max{( 0 +|pe;cpil 0 ^), 1/17) = 0.69. Using Eqs. ([T]i, ([b]), 
and (IT^ . one can arrive at 

EpigexpC) > 2.68, Ecigexpf) > 0.45, 

C(gexp,)> 0.92, yV(p„-p,)> 10.73, 

for this real experimental state. 

The second example is also from a real experiment in. 
In Ref. the authors experimentally realized a spe¬ 

cial three-photon pure state: | 0 j) = -^| 000 ) + -^| 110 ) H- 
||111) and a four-photon Dicke state with two exciations 
IDj) = -^(10011) -t- 10101) -I- 10110) -F 11001) -F 11010) + 
11100)). The square of fidelities for the special state and 
the Dicke state are measured to be - 0.9821 

and (D^lpexpsIT)^) = 0.9780, respectively. Eor simplicity, we 
only consider the entanglement under A\BC bipartition for 
the three-photon special state and AB\CD bipartition for the 
four-photon Dicke state. Other bipartition entanglement can 
also be calculated analytically based on our theorems. There¬ 
fore, our parameters for these experimental states should be 
A 2 = max{{tfs\gexp.,\ffs)/isim), 1/m} = 0.6547 withm = 2 and 
Si = 3/4, and A 3 = max{{D^^\gexp,\D^^ / (s\m), l/m) = 0.3667 
with m - 4 and si = 2/3. Using Eqs. Q, ([9]l, and ([T2l i. 
one can arrive at 


Epigexpi) > 0.1661, 
Cigexp^) > 0.3094, 
Ccigexpf) > 0.3094, 
Epigexp,) > 0.1437, 
C(gexp,) > 0.1905, 


Ecigexpf) > 0.0245, 
N(gexp,) > 0.3094, 

Ecigexp,) > 0.0160, 
Migexp,) > 0.4668, 


for these real experimental states. 

In the first example shown above, Tonolini et al. used 
compressive sensing to reduce the number of measurements, 
and 2506 measurements are needed for state reconstruction. 
Actually, the whole state reconstruction is not necessary for 
our lower bounds. We only need to perform local mea¬ 
surements (or one global measurement) to obtain the ex¬ 
pectation value of |0){0| for d x d dimensional states. If 
the global measurement |0)(0| is not easy to be performed 
in experiments, one can use local measurements instead of 
it. We recall a complete set of local orthogonal observables 
(LOOs) introduced in Ref. 1^ . i.e., {Gf} — {|/)((|, (|m)(n| + 
|«)(m|)/ V 2 , (|m)(«| - |«)(m|)/(! V 2 )), where 1 < I < d and 
1 < m < n < d. The LOOs satisfy the orthogonal con¬ 
dition Tt(Gi;Gic') = dkk', and the complete condition A - 
Yjk T'r{GkA)Gk for an arbitrary observable A. Thus, |0)(0| = 

ifLmAGkVl ® VBGlvl\(l>)VAGkVl ® VbGIvI, where G[ 
means transpose of G^ which is another complete set of LOOs, 
and {<0 |VaGj:U] ® VbGJU^I^)) = {s;, x/s„s„„ xfsfsf,}. There¬ 
fore, only d^ local measurements are required which is much 
less than the number of measurements d'^ required by tomog¬ 
raphy. Eor the first example d - 17, we only need 289 local 
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measurements, much less than 2506 measurements needed by 
the compressive sensing. When |0) = \^^), the expression of 
has been widely known from Ref. S. Moreover, 
Ref. 11^ proposed methods to get a lower bound of the maxi¬ 
mally entangled fraction with a rather smaller number of local 
measurements. 


VIII. RELATIONS BETWEEN ENTANGLEMENT 
MEASURES 


fore. 


ecio) > 




> 1-1 


> 1 


y ■> ' 

yjl -c(g)'^. 


holds, and since c(g), edg) e [0,1] we get Eq. (O. 


There have been comparative studies of entanglement mea¬ 
sures. Horodecki et al. introduced axiomatic approach for en¬ 
tanglement measures i^ . Eltschka ef a/, proposed inequali¬ 
ties between the concurrence and CREN for any m®n{m < n) 
quantum state g IH], N(g) > C{g) > yj2l[m{m- l)]A/’(p). 
We shall study the relations between concurrence, geometric 
measure of entanglement and CREN. 

Theorem 6. — Eor any m ® n {m < n) quantum state g, its 
concurrence C{g), geometric measure of entanglement Ecig), 
and convex-roof extended negativity Nig) satisfy 

Nig) < myil - Edg)) - 1, 

(mil), 

\ m J 
cigf H- (1 - edg)f < 1, 

where yix) - [^/x + yjim - 1)(1 - x)]^/m with x 6 [ 1 /m, 1 ] 
and the integer m > 2, cig) - V'«/[2(m- l)]C(p), edo) - 
mEdg)lim- 1). 

Proof. — Let {pj, \ fj)} be the optimal ensemble for Edg) 

and be the Schmidt coefficients in decreasing order for 
Thus, 


(16) 

(17) 

(18) 


m I _ 2 

Me) < = 

J J '=1 

- Yj - 1 

j 

< my{Y^pjp{)-I ^myil - Edg))- 


since yix) is a concave function and we have used the fact 
that 2 ; - -fpf + -fif < -fpf + sjim - l)i'Zi>iPi) = 

^Jpf + .y/(m2n3(T2r^. 

In order to prove Eq. ([TTl i. the properties of yix) are needed: 
yix) is a decreasing function, i.e., yix) > yiy) if x < y, and 
yiyd)) = X. Therefore, from Eq. (ITbl l we have y[Nig)lm + 
1/m] > 1 - Edg), which is equivalent to Eq. (O. 

It is worth noticing that 2,//^ = p\ + ^ + 

iIji>iPi)^lim-l) = \-2px+mjl]lim- 1 ) where/ii = 1 -pi. 
Thus, m(l - 2; pf)/itn- 1) < 2m/ii/(m- l)-m^/Ij/(m- 1)^ = 
1 - [1-m/Ii/(m-1)]^ holds. Together with pi < im-l)lm,we 

have 1 - mpilim - 1) < yjl -m(l - I], p^)/(m- 1). There- 


IX. DISCUSSION AND CONCLUSION 

Actually, the method we used to get lower bounds of EOE, 
GME, concurrence, CREN and G-concurrence can be gen¬ 
eralized to arbitrary bipartite convex-roof entanglement mea¬ 
sures. Suppose that the entanglement measure on an m ® n 
im < n) pure state |i/f) is £(|i/f)) = /(/(). where ft is the 
Schmidt vector. We first get the minimal admissible /(/?) 
for a given A = iX'JU i.e., E(T) = min^{/(p)|d = 

(2;=i ■\/pi)^ Thus, for a general m ® n im < n) state g, 
the lower bound of this entanglement measure Eig) is given 
by Eig) > co[E(A)] with A = max{(0|p|0)/(sim), 1/m), if the 
final function co[E(/l)] is a monotonously increasing convex 
function with respect to A. 

Our paper actually provided lower bounds on a variety of 
bipartite entanglement measures from a well-known entangle¬ 
ment witness OM: 


w^sii-\4>m, ( 19 ) 

where If) = 1 / 4 ® V^/Iu) is an arbitrary pure entan¬ 

gled state in m ® n system (with {-\/sJ) being its Schmidt co¬ 
efficients in decreasing order). This entanglement witness has 
been measured in many experiments, such as Refs, llb^ - l^ . 
Using this entanglement witness (where \<p) can be an arbitrary 
pure entangled state rather than only a maximally entangled 
state), one can obtain lower bounds on a variety of bipartite 
entanglement measures based on our method. 

In conclusion, we present observable lower bounds for sev¬ 
eral entanglement measures defined by convex roof, which in¬ 
clude EOE, GME, concurrence, CREN, and G-concurrence. 
The lower bounds estimate these entanglement measures for 
arbitrary finite-dimensional bipartite states. Moreover, these 
lower bounds can be easily obtained from the expectation 
value of a single observable. Based on our results, we 
present some examples using real experimental measurement 
data. In principle, our method can be used for arbitrary 
finite-dimensional bipartite convex-roof entanglement mea¬ 
sures. Last but not least, the relations between entanglement 
measures are studied. 
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APPENDIX 

Here we provide some details of calculations to get the ex¬ 
pressions of R(A), co[R(/l)], Q(A), co[2(/i)], P(A), co[f’(/l)], 
K(A) and co[fir(/i)]. The main idea is to calculate lower bounds 
of entanglement measures as a function of A for pure states, 
and then to extend the bounds to mixed states by convex hull. 
The details of proving the inequality Yjj PjA-' > A have also 
been presented. 



FIG. 1: (Color online) H(jl) versus A. 50000 red dots represent ran¬ 
domly generated states with m = 4. The lower boundary is a smooth 
blue curve, which corresponds to R(/l). 


A. Calculation of R(/l) and co[R(/l)] 


In the following, we shall seek the minimal H{(j) — 
- mil Pi log 2 /t/ for a given A = (2'" j We use R{A) 

to denote the minimal H(p) for a given A, i.e.. 


R(A) = min 

4 




(SI) 


As shown in Ref. Cl, the minimal H(/T) versus A corresponds 
to fl in the form 




l-t 
m — 1 ’ 


1 -t ] 

1 

-> for f € 

-,1 

m — 1J 

im J 


(S2) 


with m -1 copies of (1 -f)/ (m- 1) and one copy of t. Therefore, 
the minimal H{jj) and corresponding A are 

H(t) = -flogjf-)!-f)log 2 ^^— j, (S3) 

m — 1 

4(0 = i(Vf+ (S4) 

In order to show the minimal //(/O versus A, we need the in¬ 
verse function of A(t). After some algebra, one can arrive at 


B. Calculation of Q(A) and co[g(/l)] 


We first seek the minimal G{p) = 1 - for a given A - 
(2/ii where - max{yu,). We use Q(A) to denote 

the minimal G(fi) for a given A, i.e.. 


Q(A) = min 
4 




(S6) 


It is interesting that, similar to the minimal Gift) versus 
A corresponds to fl in the form 






1 -f 
OT — 1 ’ 


1 - f t 

1 

-i for f e 

-,1 

m — 1J 

im J 


(S7) 


withm-1 copies of (l-0/(»J-l) and one copy off. Therefore, 
the minimal G{p) and corresponding A are 


G(f) = 1-f, (S8) 

4(f) = l(Vf+ (S9) 

where i^max - t, since f > (1 - f)/(m - 1) with f 6 [l/m, 1]. In 
order to show the minimal G(fT) versus 4, we need the inverse 
function of A(t). After some algebra, one can see that 


f(4) = -(VI-H xJ{m-l)il-A)f, (S5) 

with 4 e [l/m, 1]. Substituting Fq. (IS5b into Fq. (IS3I) . we can 
get the expression for R{A). 

We can simulate the lower boundary of the region in H{fi) 
versus 4 plane. 50000 dots for randomly generated states with 
m = 4 are displayed in Fig. [1] The lower boundary corre¬ 
sponds to R(A). It is defined that co(g) is the largest convex 
function that is bounded above by the given function g. From 
the expression of R(A), the explicit expression of co[R(A)] is 
obtained in Refs. which is Fq. (2) shown in the 


f(4) = -(V4-H J(m-l)(l-A)f, (SIO) 

with 4 6 [1/m, 1]. Substituting Fq. (ISlOl l into Fq. (IS8I) . we 
can get the expression for Q(A), i.e., 

Q(A) = 1 - VI + - 1)(1 - 4)f. (Sll) 

From Fq. (ISl lb . we can find that 

dQ(A) _ ( “ ^)( ^ V(fn- 1)(1 -4)) ^ ^ 

~d4 


main text. 


m 
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FIG. 2: (Color online) G(jl) versus A. 50000 red dots represent 
randomly generated states with m = 4. The lower boundary is a 
smooth blue curve, which corresponds to both Q(A) and co[2(/l)], 
since co[2(d)] = Q(A). 


when A e [l/m, 1]. Thus, Q{A) is a monotonously increasing 
function. Moreover, from the definition of co(), one can see 
that co[2(/l)] = Q(A), since Q(A) is a convex function. In 
order to prove this, we only need to show that ^Q{A)ldA^ > 0. 
From Eq. dSl 11 1, one can get 


d^Q(A) V(m-l)(l-T) 

dT2 2 ot(1-T)2t3/2 - ’ 


(S12) 


since m is an integer (m > 2) and A 6 [1/m, 1]. Therefore, 
co[2(T)] = Q(A), which is a monotonously increasing convex 
function. 

We can simulate the lower boundary of the region in G(jl) 
versus A plane. In Fig. |2] 50000 dots for randomly generated 
states with m = 4 are displayed. The lower boundary corre¬ 
sponds to Q(A), which coincides with co[2(/l)]. 


withm-l copies of (I-f)/(m-l) and one copy off. Therefore, 
the minimal L(Jl) and corresponding A are 

£.(,) = (SI 5) 

V m — 1 

A(t) = -(V^+ V(1 -0(fn- l)f. (S16) 

In order to show the minimal L{p[) versus A, we need the in¬ 
verse function of /1(f). After some algebra, one can arrive at 


t{A) = ^(Vd-t- y/im-l)(l-A)f, (S17) 

with A 6 [1/m, 1]. Substituting Eq. (IS17b into Eq. (IS151 l. we 
can get the expression for P{A), i.e.. 


PU) = 


2(1 — t){mt + m — 2) 


m — 1 


f = 1(VI-H y/im-m-A)f. 


From Eqs. (ISlSb and (IS 19b . we can find that 


(SI 8) 
(SI 9) 


dE(/l) _ dP(f) df(d) 
dA df dA 

where 


(S20) 


dP(0 

df 

df(d) 

dd 


V2(l-mf) ^ 

^(m — 1)(1 — f)(m + mt — 2) 
1/1 1 - m \ 

V(m- l)(l-i5j 
x( Vd -I- V('«-1)(1 -d)) < 0, 


(S21) 


(S22) 


since m > 2 and f, d e [1/m, 1]. Therefore, 


dP(d) 

dd 


> 0 , 


(S23) 


which means P(A) is a monotonously increasing function. 

Furthermore, in order to show P(A) is a concave function, 
we need to prove that d^P(d)/dd^ < 0. One can use 


C. Calculation of P(A) and co[P(d)] 


We first seek the minimal L(/I) = ■^2(1 - Yj'ILi pj) for ^ 

given d = iY/iU 1^- We use P(A) to denote the minimal 
L(/j) for a given d, i.e.. 


P(A) = min < L(/J) 


d = 


1=1 '' 


(SI 3) 


It is interesting that, similar to N(/I) and G(/f), the minimal 
L(p) versus d corresponds to (1 in the form 


— 

1 m — 1 


1 -f 
’ m — 1 


for f e 


-,1 

m 


(S14) 


d^P(A) _ d2p(f)/df(d)\^ dP{t)d^t{A) 
dd^ df^ \ dd / '*' df dd^ 


(S24) 


and 

d^Pit) 

dfi 

dh{A) 

dd 2 

Thus, 


-m-2-0/2^ 

HI — t)(m + mt — 2)' 


V(m- 1)(1 -d) 
2 m(l - d)2d3/2 ■ 


(525) 

(526) 


d^P{A) 

dd2 

b‘^{m — 2)[(m — \)a — b]^f 
V2m2(m — l)^^2(i _ d)2d2/2(l — t)^G(m + mt — 2)^/2 





















































8 



FIG. 3: (Color online) L(jl) versus A. 50000 red dots represent ran¬ 
domly generated states with m = 4. The lower boundary is a smooth 
blue curve which corresponds to P(A), and the dashed black line cor¬ 
responds to co[f’(/l)]. 



FIG. 4: (Color online) S (jT) versus A. 50000 red dots represent ran¬ 
domly generated states with m = 4. The lower boundary is a smooth 
blue curve which corresponds to K(A), and the dashed black line cor¬ 
responds to co[A'(/i)]. 


where 


a Vt, 

b -\J(m — 1)(1 — x), 
a 

X 

f —3{m — l)x^ + 2(3 — m)x + 


2m — 3 


m — 1 


with X 6 [l/(m - 1), -Hoo). It is easy to see that 

f<0, 

when X e [l/{m - 1), -i-oo). Thus, 


dT2 


< 0, 


(527) 

(528) 

(529) 

(530) 

(531) 

(532) 


when OT >2,t,Ae [l/m, 1], and x 6 [1/(ot- 1),-hoo). There¬ 
fore, the convex hull of P(A) will be 


co[P(T)] = 


m — 1 m 


(S33) 


which is a straight line from (1/m, 0) to (1, V2(m - l)/m). 

We can simulate the lower boundary of the region in L(/() 
versus A plane. 50000 dots for randomly generated states with 
m - 4 are displayed in Fig. [3 The lower boundary cor¬ 
responds to P{A), and the dashed black line corresponds to 
co[F(T)]. 


D. Calculation of K(A) and co[A'(/l)] 


As shown in Ref. ifs^ . the minimal S (jl) versus A corresponds 
to /? in the form 




1 -f 
’ m — 1 ’ 


1 -f 
m — 1 


for t e 



(S35) 


withm-1 copies of(l-f)/(m-l)and one copy off. Therefore, 
the minimal S (ft) and corresponding A are 

/ (1 - f)"’-^ 

A(t) = -(Vf+ (S37) 

In order to show the minimal S (/t) versus A, we need the in¬ 
verse function of A(t). After some algebra, one can arrive at 

t(A) = i(VI- ^(m-l)(l-A)f, (S38) 

with/l 6 [(m-l)/m, 1]. Substituting Eq. ( IS38I ) into Eg. ( IS36I ), 

we can get the expression for K{A), i.e., 

K(A) = m[f(l(S39) 

f = -(^^A- ^J(m-l)(l-A)f. (S40) 

Erom Ref. JUI, one can see that K{A) is a monotonously 
increasing concave function in [(m - l)/m, 1]. Therefore, the 
convex hull of K(A) will be 


We first seek the minimal S {ft) = m(n™ i for a given 
^ = {YjT=i yfJfif'I^- We use K{A) to denote the minimal S(ft) 
for a given A, i.e.. 


K{A) - imn \S(ji) 



(S34) 


co[/:(T)] = max{l-m(l-T),0). (S41) 

We can simulate the lower boundary of the region in S (ft) 
versus A plane. 50000 dots for randomly generated states with 
m - 4 are displayed in Eig. H] The lower boundary cor¬ 
responds to K(A), and the dashed black line corresponds to 
co[/:(T)]. 
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E. Proof of the inequality 2^ pjA^ > A 

Here we give the details of proving 2; Pj'^^ ^ For any 
m®n {m < n) quantum state g, suppose that we have found an 
optimal decomposition Yjj forp to achieve the infi- 

mum of E{q), where E{g) is one kind of entanglement mea¬ 
sure defined by the convex roof. For each \4<j), we have the 
expression 

m m n 

|i/r/) = Ua®UbY^ -y/t/lu) ^ 

i=l i=l k=l 

with { being its Schmidt coefficients in decreasing order, 
where Ua (Ub) is an m x m (n x n) unitary matrix and S is 


an m X n matrix defined by = ^p\5ik. Similarly, for an 
arbitrary pure entangled given state \(p) in m ® n system, we 
have the expression 

m m n 

\^) = Va ® Fb ^ Vs;|//) ^Va®Vb ZZ Siklik), (S42) 

1=1 1=1 it=i 

where {-^1 are its Schmidt coefficients in decreasing order, 
Va (Vb) is an m x m (n x n) unitary matrix, and S is an m x n 
matrix defined by S ik - ylTidik- Therefore, the maximum 
under all possible U\ (which is an m x m unitary matrix) and 
U 2 (which is an n X n unitary matrix) is 


max 

U^,U2 


{4>\Ui®U2\iiJj){iifj\u\®ul\4>) 

siin 


max 

UuVi 


12::;v=i n,k'=i ^ikSi,,{i'\viuiUA\i){k'\viu2UB\k)\^ 


siin 


max 

Ul,U2 


< max 

UuUi 


|Tr('l/iZ'Z/j5')p 

s\m 

i:;"i 

Siin 


s\m 



s\m 


(S43) 


where Hi is an m x m unitary matrix defined by 'Z/i,-,- = {i'\V\U\UA\i), H 2 is an « x n unitary matrix defined by 'U2k'k - 
{k'\V"kU 2 UB\k), and {cr,(A)) are singular values of A in decreasing order. The first inequality holds since |TrA| < ZiCriCA) for 
arbitrary matrix A, and the second inequality holds since the following theorem S: 

Let A (nx p matrix) and B (pxm matrix) be given, let q - mm{n, p, m], and denote the ordered singular values of A, B, and 
AB by o"!(A) ^ ^ ^ 0, o~i(^B) ^ ^ *^min{/j,m}(^) — Oj and o~i(AB) ^ ^ — 0. Then 


^ cr,(AB) < ^ cr,(A)cr,(B). 

i=l i=l 


(S44) 


Therefore, 


A-' = 


(zzi (i:;=i 


> 


A 2 


m s\ni 

where the first inequality holds since si = max{s,}. Thus, 


> max 

Ul,U2 




{cl2\U, ® U2\ilfj){ilfj\u\ ® t/>) 

S\m 


{^\Uy ® U 2 Y 2 jPjWj){'kj\U\ ® t/210) {^\q\cI)) 

max-> -. 

Ui,U2 Sim sim 


Together with A^ - {YIZi ^ ^Im, one can obtain Pyd-' > vaax{{^\Q\(p)l{sim), l/m) = A. 


(S45) 


(S46) 
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